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Particle velocity distributions measured in the weakly collisional solar wind are frequently found 
to be non-Maxwellian, but how these non-Maxwellian distributions impact the physics of plasma 
turbulence in the solar wind remains unanswered. Using numerical solutions of the linear dispersion 
relation for a collisionless plasma with a bi-Maxwellian proton velocity distribution, we present a 
unified framework for the four proton temperature anisotropy instabilities, identifying the associ¬ 
ated stable eigenmodes, highlighting the unstable region of wavevector space, and presenting the 
properties of the growing eigenfunctions. Based on physical intuition gained from this framework, 
we address how the proton temperature anisotropy impacts the nonlinear dynamics of the Alfvenic 
fluctuations underlying the dominant cascade of energy from large to small scales and how the 
fluctuations driven by proton temperature anisotropy instabilities interact nonlinearly with each 
other and with the fluctuations of the large-scale cascade. We find that the nonlinear dynamics of 
the large-scale cascade is insensitive to the proton temperature anisotropy, and that the instability- 
driven fluctuations are unlikely to cause significant nonlinear evolution of either the instability-driven 
fluctuations or the turbulent fluctuations of the large-scale cascade. 


I. INTRODUCTION 

The near-Earth solar wind is a dynamic plasma 
environment supporting broadband turbulent spec¬ 
tra of plasma and electromagnetic field fluctuations, 
providing a uniquely accessible venue for the study 
of the fundamental physics of astrophysical plasma 
turbulence. Direct spacecraft measurements show 
that the particle velocity distributions of ions and 
electrons in the solar wind commonly deviate from 
the isotropic Maxwellian velocity distributions char¬ 
acteristic of a plasma in local thermodynamic equi¬ 
librium, a result that is not unexpected considering 
the weak collisionality of the solar wind plasma. But 
how the non-Maxwellian nature of the plasma par¬ 
ticle distribution functions impacts the physics of 
plasma turbulence in the solar wind remains unan¬ 
swered. Unraveling the nature of turbulence in the 
solar wind is a grand challenge problem in helio¬ 
physics because turbulence significantly impacts the 
transport of energy from large-scale motions to suffi¬ 
ciently small scales at which that energy is efficiently 
converted to plasma heat or some other non-thermal 
form of particle energization. 

Turbulence in the solar wind is dominated by en¬ 
ergy injected into the turbulent cascade at large 
scales through nonlinear interactions among plasma 
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motions at those large scales. The turbulent cascade 
from those large scales down to characteristic ion ki¬ 
netic scales, denoted the inertial range of the tur¬ 
bulence, occurs without significant dissipation and 
is mediated by nonlinear interactions between in¬ 
compressible Alfvenic fluctuations with wavevectors 
that become increasingly perpendicular (fcj_ ^ fc||, 
where perpendicular and parallel are defined relative 
to the local mean magnetic field) with diminishing 
length scale. In addition, this turbulent cascade of 
energy from large scales through the inertial range 
also includes a small admixture of compressible fluc¬ 
tuations, which appear to have similarly anisotropic 
wavevectors with k±_ ^ ^||'^*^ Separate from this 
anisotropic cascade of energy from large scales (re¬ 
ferred to as the large-scale cascade in this paper), ki¬ 
netic instabilities driven by non-Maxwellian velocity 
distributions may also inject energy directly into the 
turbulence typically at small scales near the charac¬ 
teristic ion or electron kinetic scales. 

As the first step of an in-depth study of the im¬ 
pact of non-Maxwellian distribution functions on 
plasma turbulence in the solar wind, we focus here 
specifically on the effect of a bi-Maxwellian pro¬ 
ton temperature distribution. We investigate two 
open questions: (1) How does the proton temper¬ 
ature anisotropy impact the nonlinear dynamics of 
the Alfvenic fluctuations underlying the large-scale 
cascade? (2) What is the nature of the modes gen¬ 
erated by kinetic instabilities driven by the proton 
temperature anisotropy, and how do they contribute 
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to the fluctuations measured in the solar wind? 

In Section ini we discuss the physics of the proton 
temperature anisotropy and the resulting kinetic in¬ 
stabilities, and we describe the numerical approach 
used to solve for the linear kinetic plasma behavior. 
A unified framework for understanding the four pro¬ 
ton temperature anisotropy instabilities is presented 
in Section Ell The physical properties described in 
this framework are employed to discuss the effect 
of the proton temperature anisotropy on the large- 
scale turbulent cascade in Section lTVl and to examine 
how fluctuations generated by the kinetic instabili¬ 
ties interact nonlinearly with each other and with 
the turbulent fluctuations of the large-scale cascade 
in Section E Our findings are summarized in Sec¬ 
tion ED and a brief technical note on the identifica¬ 
tion of normal modes of the Vlasov-Maxwell system 
appears in Appendix El 


II. THE PHYSICS OF PROTON 
TEMPERATURE ANISOTROPY 

Compressions or expansions of a magnetized 
plasma result in different rates of change of charged 
particle velocities in the directions perpendicular 
and parallel to the local magnetic field. Collisions 
act to isotropize the velocity distributions in differ¬ 
ent directions, ultimately driving the velocities to 
a Maxwellian distribution, corresponding to local 
thermodynamic equilibrium with a single isotropic 
temperature, a state of maximum entropy from 
which no free energy can be extracted. But, in the 
solar wind, collisions are often insufflcient to yield 
isotropy between parallel and perpendicular temper¬ 
atures, representing a potential source of free energy 
in the particle velocity distributions. If the plasma 
motions lead to temperatures that are sufficiently 
anisotropic, then kinetic temperature anisotropy in¬ 
stabilities can tap this source of free energy to drive 
electromagnetic fluctuations in the plasma that ulti¬ 
mately serve to reduce the temperature anisotropy. 

Note that other types of non-Maxwellian velocity 
distributions, such as the presence of a beam com¬ 
ponent or a relative drift between plasma species, 
also contain free energy that may drive kinetic in¬ 
stabilities. But, to build a foundation upon which 
to understand the effect of non-Maxwellian velocity 
distributions on plasma turbulence, we begin with 
the idealized case of a bi-Maxwellian proton temper¬ 
ature distribution, in which the free energy content 
is characterized by the single parameter Tlp/Tjip. 

Spacecraft measurements in the near-Earth solar 


wind demonstrate a wide spread of values of the 
proton temperature anisotropy lAp/Tup, often fill¬ 
ing the entire range between the marginal stability 
boundaries of kinetic proton temperature anisotropy 
instabilities^^— The physical cause for this spread 
of Tj_p/T||p values remains to be definitively de¬ 
termined. For double adiabatic evolution of the 
temperatures in a magnetized plasma?^ the spher¬ 
ical expansion of the solar wind in the inner helio¬ 
sphere with a typical Parker spiral magnetic field 
leads to Tj_p/r||p < 1— On the other hand, 
mechanisms proposed to yield Tj_p/r||p > 1 include 
proton cyclotron heatingjiS shocks^ compression 
between solar wind streams)^, compressional slow 
wave modeSfi^ or double adiabatic expansion with a 
transverse magnetic field— A discussion of the evo¬ 
lution of T±p/T||p with heliocentric radius can be 
found in Matteini et. al. 2007—, and a review of the 
many potential mechanisms governing TLp/T||p and 
their associated time scales in the solar wind is pre¬ 
sented in TenBarge et. al. 2015^^. Here we will not 
further address the causes of the proton temperature 
anisotropy, but merely accept it as an observational 
fact and explore the consequences for the turbulence. 

For a plasma with a bi-Maxwellian proton distri¬ 
bution and plasma parameters relevant to the solar 
wind, there exist at least four potential electromag¬ 
netic proton temperature anisotropy instabilities: 

( 1 ) the parallel (or whistler) firehose instability 

(2) the Alfven (or oblique) firehose instabilityj^i (3) 
the mirror instability^^^— and (4) the proton cy¬ 
clotron instability— The first two of these instabil¬ 
ities occur for T^p/T||p < 1 and /?||p > 1, while the 
latter two instabilities arise for Ttp/Typ > 1 and all 
/3||p, where /3||p is the ratio of parallel proton thermal 
pressure to magnetic pressure. When the plasma 
exceeds a threshold value of the proton temperature 
anisotropy, these instabilities can tap the free energy 
associated with the anisotropic proton temperature, 
driving electromagnetic fluctuations and ultimately 
reducing the temperature anisotropy, thereby mov¬ 
ing the plasma back toward a state of marginal 
stability.—^— 

It is important to emphasize that each of these 
instabilities is necessarily associated with a normal 
wave mode, or eigenmode, of the kinetic plasma. 
Each solution of the linear kinetic dispersion relation 
yields a complex eigenfrequency cj + ij, where the 
sign of the imaginary component indicates whether 
the mode is growing or damped. For a uniform 
plasma with an isotropic Maxwellian distribution, 
all eigenmodes are damped. But if the proton tem¬ 
perature is sufficiently anisotropic, unstable modes 
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may grow, utilizing the free energy in the velocity 
distribution. This leads to the injection of energy 
into electromagnetic fluctuations in the solar wind 
plasma at scales where the instability growth rate 
is positive. Here we aim to make clear the connec¬ 
tion between each instability and its associated wave 
mode, to elucidate the properties of these instability- 
driven modes, and to discuss the resulting contribu¬ 
tion to the fluctuations measured in the turbulent 
solar wind. 


A. Numerical Solution 

To explore the properties of a plasma with a bi- 
Maxwellian proton velocity distribution, we employ 
the numerical Vlasov-Maxwell linear dispersion re¬ 
lation solver PLUME (Plasma in a Linear Uniform 
Magnetized Environment). PLUME extends the 
solver described in Quataert 199821 by allowing for a 
bi-Maxwellian equilibrium temperature distribution 
for both electrons and an arbitrary number of ion 
species. We have benchmarked our results against 
the widely used linear dispersion relation solvers 
WHAMP— and NHDSH and found agreement. 

The general linear dispersion relation for a fully 
ionized, proton-electron plasma with bi-Maxwellian 
particle distributions can be expressed as 

_ (l 1 ^ a ^IIp 

WVM — P||p, I rp J rp I )■ 

UIp U|e U|e C 

( 1 ) 

The equilibrium magnetic field is Bq = and 

we solve for the eigenfrequencies and eigenfunctions 
of a plane wave with wavevector k = k±i -|- fc||Z. 
The proton gyroradius and inertial lengths are de¬ 
fined as Pp = Vt±p/i^p and dp = cjujpp, the per¬ 
pendicular and parallel temperatures (expressed in 
units of energy) for species s are T±s and Tp, the 
proton thermal velocities parallel and perpendicu¬ 
lar to Bo are = y^2r(||^_L)p/mp, and the 

speed of light is c. The parallel proton plasma 
beta, or ratio of parallel thermal to magnetic pres¬ 
sure, is /3||p = 87rnT||p/i?o, and the proton gyrofre- 
quency and proton plasma frequency are given by 
Hp = qB^/irUpc) and ujpp = yj^impq^/nip. 

In order to focus on the four proton temper¬ 
ature anisotropy instabilities, we consider a non- 
relativistic plasma, Vt\\p/c <C 1, in which the elec¬ 
trons have a Maxwellian equilibrium distribution, 
T^e/T\\e = 1; nnd the protons and electrons have 
equal parallel temperatures, T||p/T||e = 1. The 
proton distribution is bi-Maxwellian, with distinct 
perpendicular and parallel temperatures, Tj^p and 


T\\p. Under these assumptions, linear solutions of 
the Vlasov-Maxwell dispersion relation are depen¬ 
dent on only four parameters, 

T 

wvM = a;(fc_LPp,fc||dp,/3||p, 7^). (2) 

^ lb 

We avoid more complex distributions—such as 
a Kappa distribution)21“i2£ distributions composed 
of several overlapping Maxwellians,— >22 or multiple 
ion species—because the non-Maxwellian nature of 
the bi-Maxwellian solution depends only on the sin¬ 
gle parameter Tup/Tup, providing physical insight 
into the fundamental physics of these instabilities 
without the complications introduced from more 
complicated non-Maxwellian distributions. Future 
work concerning the alteration to the temperature 
anisotropy instabilities from other non-Maxwellian 
distributions^^— additional ion species^ and rel¬ 
ative drifts2i2£i^ will prove useful in extending the 
results of this work to achieve a more complete char¬ 
acterization of the behavior of turbulent fluctuations 
in the solar wind. 

One of our main aims is to establish a clear 
connection between the four proton temperature 
anisotropy instabilities and the low-frequency nor¬ 
mal modes of the Vlasov-Maxwell system, specifi¬ 
cally the kinetic counterparts of the MHD Alfven, 
fast, slow, and entropy modes.— Each of these four 
normal modes is associated with a distinct dispersion 
surface . A dispersion surface is a map formed 
by the solution to the linear dispersion relation (for 
example, the real component of the complex eigen- 
frequency) over wavevector space, in this case the 
(fc_L, fc||) plane. Different regions in wavevector space 
on a single dispersion surface correspond to different 
commonly known wave modes that may have dis¬ 
tinct properties. 

To avoid confusion, we specify here each of these 
four dispersion surfaces (associated with the Alfven, 
fast, slow, and entropy modes) and identify the com¬ 
monly known wave modes represented by different 
regions on each surface. The Alfven dispersion sur¬ 
face includes the MHD Alfven wave at k±pp <C 1 
and fcydp <C 1, the kinetic Alfven wave at k^pp>l 
and k\\dp <C 1, and the proton cyclotron wave at 
^ 1- The fast dispersion surface includes the 
fast magnetosonic wave at k±pp <C 1 and k\\dp 1, 
the ion Bernstein wave at k±pp > 1 and fc[|dp <C 1, 
and the whistler wave at k\\dp > 1. The associations 
between these commonly known wave modes and the 
Alfven and fast dispersion surfaces are illustrated 
in Figure 1 of Howes et al. 2014— The slow and 
entropy dispersion surfaces include the slow magne¬ 
tosonic wave and non-propagating entropy mode at 
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k±Pp 1 and fc||(ip -C 1; neither of these two dis¬ 
persion surfaces has commonly used names at the 
small, kinetic scales. In the remainder of this pa¬ 
per, we choose to identify the appropriate dispersion 
surface by the associated large-scale, low-frequency 
mode: the Alfven wave, fast wave, slow wave, or en¬ 
tropy mode. A detailed discussion regarding con¬ 
necting MHD and Vlasov-Maxwell normal modes 
can be found in Krauss-Varban et al. 1994^ and 
Klein et al. 2012)^ 

We employ the following rigorous procedure 
to connect each of the four proton temperature 
anisotropy instabilities to one of the four low- 
frequency, normal modes of the Vlasov-Maxwell 
system. Consider the case that our numerical 
solver finds an unstable mode (a solution for the 
complex eigenfrequency that has a positive imag¬ 
inary component) for a given set of the four pa¬ 
rameters, V* = [{kj_pp)*,mdp)*,l3p,{Tj_p/Ty)*]. 
First, we begin with the solutions for each 
of the four normal wave modes for the initial 
parameters fc||dp,/3||p, 7j_p/r||p) = 

(10“^, 10“^, 1,1). Next, we perform a nearly con¬ 
tinuous variation of each of these four parameters 
until we reach the desired point in parameter space 
V*. The nearly continuous variation for a given 
parameter involves many small increments of that 
parameter; after each small increment, the roots 
of the dispersion relation are recalculated, using 
solutions from the previous parameters as initial 
guesses for the new solutions. Finally, the unsta¬ 
ble mode is identified as one of the four normal¬ 
mode solutions at V*. This procedure allows for 
a smooth connection between the well-established, 
isotropic large-scale normal modes and the smaller 
scale, anisotropic modes and instabilities. 

Identifying these modes by such a continuous vari¬ 
ation can be made difficult by the presence of excep¬ 
tional points on the solution surface of the Vlasov- 
Maxwell dispersion relation^^^— a technical point 
discussed in Appendix A. For example, at j5\\p > I, 
the whistler wave is actually the extension of the 
slow modej^i^ rather than the fast mode, due to a 
large-scale mode exchange described in Appendix A. 


B. The Properties of the Proton 
Temperature Anisotropy Instabilities 

As discussed in Section III Al the linear disper¬ 
sion relation for a non-relativistic, fully ionized, 
proton-electron plasma with a bi-Maxwellian pro¬ 
ton distribution depends only on the four param¬ 


eters (fcj_Pp, fc||dp,/3||p, Tlp/T'iip). Therefore, the 
two plasma parameters, f3\\p and Tip/Tup, con¬ 
trol whether the plasma is unstable. In the 
event of instability, some regions of wavevector 
space {k±pp^k\\dp) will have positive growth rates, 
where observations of the electromagnetic fluctua¬ 
tions driven by the instability are generally expected 
to be dominated by modes with wavevectors having 
the maximum growth rate. The four proton tem¬ 
perature anisotropy instabilities arise in different re¬ 
gions of (/3||p, rj_p/T||p) plasma parameter space and 
lead to unstable wave growth in distinct regions of 
(fc_L/ 0 p, fcydp) wavevector space. Determination of 
these regions provides useful information concern¬ 
ing how these instabilities alter the existing large- 
scale turbulent cascade and inject energy into the 
turbulence in the solar wind, the key topics of this 
work. In Section lllll we classify the underlying fun¬ 
damental mechanisms which control each of the four 
proton temperature anisotropy instabilities, and use 
the properties illustrated in Figures [2 |3l 01 and [5] 
to examine the behavior of the four modes gener¬ 
ated by these instabilities as a function of the four 
parameters (fcj_Pp, k\\ dp, /?||p, T^_p/T\\p). 

Hellinger et al. 2006^ has compiled values for the 
marginal stability boundaries for these instabilities 
in the (/3||p, rLp/T||p) plane. The marginal stabil¬ 
ity boundary is set by calculating a;(k) for a fixed 
/3||p over all possible wavevectors k. The proton 
temperature anisotropy T±p/T\\p is then varied until 
the most unstable wavevector has a growth rate of 
■j/flp = 10“^, thus establishing the instability cri¬ 
terion. These criteria are generally well fit by an 
expression of the form 

‘ + (/3||, - M‘ 

where a, b, and /3o are unique values for each of the 
four instabilities. Values from Hellinger et al. 2006^ 
are found in Table Hand used to plot the four thresh¬ 
olds in Figure [2 Observations of /3||p and Tlp/T||p 
are approximately constrained by the mirror and 
Alfven firehose boundaries in the solar wind4“— 

For a single point in the {j3\\p,T±p/T\\p) plasma 
parameter space at which one of the proton tem¬ 
perature anisotropy instabilities arises, the unsta¬ 
ble modes will only occupy a particular region 
of {k±pp, k^^dp) wavevector space. Generally, the 
parallel firehose and proton cyclotron instabilities 
have their largest growth rates for nearly parallel 
wavevectors with k^^dp ^ 1 and k±pp <C 1, while 
the mirror and Alfven firehose are most unstable at 
oblique angles with fc||dp ~ k±pp ^ I. To illus¬ 
trate the unstable regions of {k±pp,k\\dp) wavevec- 
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FIG. 1. Marginal linear stability thresholds for the four 
proton temperature anisotropy instabilities using Equa¬ 
tion [3] and the parameters in Table H] 
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TABLE I. Instability Threshold Parameters for 
= 10“® from Helliuger et al. 2006^. 



a 

b 

do 

Proton Cyclotron 

0.43 

0.42 

-0.0004 

Parallel Firehose 

-0.47 

0.53 

0.59 

Alfven Firehose 

-1.4 

1.0 

-0.11 

Mirror 

0.77 

0.76 

-0.016 


tor space for each of the four instabilities, we present 
in Figure [5] maps of the positive growth rates 'y/flp 
for the modes associated with the four instabili¬ 
ties in (fc_LPp, fc||dp) wavevector space at fixed val¬ 
ues of /?||p and Tj_p/T||p. A discussion of the clas¬ 
sification of microscale instabilities, with k±pp ~ 1 
and/or fc||dp ~ 1, and macroscale instabilities, with 
kj_Pp <C 1 and fc||dp ^ 1, is presented in Section Hill 
Note that the relation between the proton inertial 
length and gyroradius depends on both the parallel 
plasma beta and the proton temperature anisotropy, 
Pp = dp^y /3||pTLp/T||p. 

To visualize the connection between each proton 
temperature anisotropy instability and its associated 
linear eigenmode, we present in Figure [3] solutions 
for the fast (red), Alfven (black), slow (blue), and 
entropy (gray) modes in complex frequency space 
(w/rip, y/rip) parametrized as a function of Tj_p/T[|p. 


For each panel, we chose appropriate values of /3||p, 
kj_Pp and k^^dp to allow for the development of one of 
the proton temperature anisotropy instabilities and 
vary Tj_p/T||p from unity, indicated by an open circle, 
to Tlp/Tiip = 10 (Tip/Tjip = 0.1) indicated with an 
open diamond (triangle) in the right (left) panels. 
The oblique instabilities, with k±pp ^ fc||dp ~ 1, are 
plotted in the top row and the parallel instabilities, 
with k±pp <g; k\\dp ~ 1, are plotted in the bottom 
row. 

Figure [3] shows clearly that the proton cyclotron 
instability (lower right) is connected to the proton 
cyclotron wave (black), and is therefore associated 
with the Alfven wave. The mirror instability (upper 
right) is connected to the non-propagating entropy 
mode (grey). The parallel firehose instability (lower 
left) is connected to the whistler wave (red), and is 
therefore associated with the fast wave. Finally, the 
Alfven firehose instability (upper left) is connected 
to non-propagating Alfvenic fluctuations (black). It 
should be noted that the unstable Alfven wave solu¬ 
tion (shown in black) for the mirror instability pa¬ 
rameters (upper right) is due to an extension of the 
proton cyclotron instability to oblique angles for the 
very anisotropic temperature ratio T±p/T\\p ~ 10, 
and is not related to the mirror instability. The be¬ 
havior of each of these parametrized paths, a key 
result of this paper, will be discussed thoroughly in 
Section Ell 

Figure |4] plots the roots of the Alfven, fast, 
slow, and entropy modes in complex frequency space 
{uj/flp, y/flp) for fixed values of /3||p, k\\dp, and k±pp. 
The values of these three parameters are distinct 
for each instability, and have been chosen to max¬ 
imize the growth of the associated unstable mode. 
The instabilities are ordered by row as proton cy¬ 
clotron (top), parallel firehose (second), Alfven fire¬ 
hose (third), and mirror (bottom). The proton tem¬ 
perature anisotropy is varied from isotropy (left col¬ 
umn) to a value near marginal stability (center), 
to a value of unstable growth (right). Also plotted 
are contours for the amplitude of the linear disper¬ 
sion relation |T>(a;, 7 )| (c.f. Equation 10-73 in Stix 
1992^), where solutions to the linear dispersion re¬ 
lation are found at points where |27| = 0. These 
contours, while not physically meaningful, do pro¬ 
vide a guide for the eye to where the solutions of the 
linear dispersion relation occur in complex frequency 
space. 

In Figure El we plot the complex eigenfrequen- 
cies and magnetic and density eigenfunctions for 
the four modes associated with the proton tempera¬ 
ture anisotropy instabilities as a function of T±p/T\\p 
and j3\\p at fixed points in k\\dp and k±pp, show- 
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FIG. 2. Maps in wavevector space of the growth rate 7 /IIp for the unstable modes of the proton temperature 
anisotropy instabilities for fixed values of /3||p and Txp/Tup- Both a marginally unstable (left column) and strongly 
unstable (center) case is presented for the proton cyclotron (first row), parallel firehose (second), Alfven firehose 
(third), and mirror (fourth) instabilities. A thin black line, indicating k±pp = fey dp, is included as an aid in distin¬ 
guishing between parallel and oblique unstable modes. In the right column are presented the points in (/3||p, Txp/Typ) 
parameter space from which these plots are generated as well as the associated marginal instability lines as calculated 
by Equation |3] 


ing the eigenfunction characteristics for both sta¬ 
ble and unstable values of the plasma parameters 
(/?||p) 7j-p/T'||p)- As with FigurelH a single wavevec¬ 
tor {k±pp, fc||(ip) is chosen for which unstable modes 
exist for some values of /3||p and T±p/T\\p. The 
instabilities are plotted by column as proton cy¬ 


clotron (first), parallel firehose (second), Alfven fire¬ 
hose (third), and mirror (fourth). The parameters 
7j_p/r||p and /3||p are both varied from 0.1 to 10, 
with the Txp/Tiip variation shown on the horizon¬ 
tal axis and the (3\\p variation shown in color, with 
blue indicating /3||p = 0.1, black indicating j3\\p = 1, 
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FIG. 3. Vlasov-Maxwell dispersion relation solutions parametrized by a varying T±p/T\\p for fast (red), Alfven 
(black), slow (blue), and entropy (gray) modes in complex frequency space (cu/flp, y/flp). Solutions in each of the 
panels represent choices of plasma parameters which will lead to the four proton temperature anisotropy instabilities: 
the Alfven firehose (top left), the mirror (top right), the parallel firehose (bottom left), and the proton cyclotron 
(bottom right). Tj^p/Typ <!(>!) modes are shown in the left (right) column and the oblique (parallel) modes are 
shown in the top (bottom) row. 


and red indicating /3||p = 10. By row from the 
top, we plot uj/flp (first), | 7 /f^p| (second), \SBx/SB\ 
(third), \SBy/6B\ (fourth), \6Bz/SB\ (hfth), and 
l^rip/nol (sixth). In the second row, the unstable 
(stable) inodes are indicated by solid (dashed) lines 
for y/fip > 0 (< 0). 


III. TAXONOMY OF PROTON 
TEMPERATURE ANISOTROPY 
INSTABILITIES 

The classification of plasma instabilities is an im¬ 
portant but often subtle matter. A variety of dif¬ 
ferent schemes have been applied, contrasting fluid 
vs. kinetic instabilities, macroinstabilities vs. mi¬ 
croinstabilities, or configuration-space vs. velocity- 
space instabilities. These different schemes classify 
the instabilities according to different criteria, such 
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FIG. 4. Plots of the linear Vlasov-Maxwell dispersion relation 7 /flp)| illustrating the transition of modes 

associated with the proton temperatnre anisotropy instabilities from stable (left column) to marginally stable (center) 
to unstable (right). Contours of constant value are shown as black lines, with eigenfrequency solutions for the Alfven, 
fast, entropy, and slow modes corresponding to |I>(a;/f2p,y/flp)] = 0 shown as red points. The plasma parameters 
used for each contour map have been chosen to generate the proton cyclotron (first row), parallel firehose (second), 
Alfven firehose (third), and mirror (fourth) instabilities. The monotonic change of rLp/Tup is the only variation of 
parameters made in each row. Discussions of the behavior of the instabilities are found in subsections [IIIlA-D. 


as the spatial scales at which the instability op¬ 
erates or the nature of the underlying mechanism 
driving the instability. For instance, Treumann & 
Baumjohann 1997— categorize instabilities as either 
macroinstabilities or microinstabilities depending on 


the spatial scales associated with the instability. 
Macroinstabilities occur at scales much larger than 
the particle kinetic scales (ps or ds), while microin¬ 
stabilities arise near these kinetic scales. Krall & 
Trivelpiece 1973^ alternatively define configuration- 




















































































































































































































































space instabilities to be those which are associated 
with the departure of macroscopic quantities from 
thermodynamic equilibrium, while velocity-space in¬ 
stabilities are those which depend on departures 
from an isotropic Maxwellian velocity distribution 
function. The configuration-space vs. velocity-space 
distinction is also often referred to as a fluid vs. ki¬ 
netic distinction. A detailed discussion about the 
relation between a plasma’s macroscopic quantities 
and associated microscopic instabilities can be found 
in Schekochihin et al. 2010^. 

It must be noted that the Treumann and Krall 
definitions are not synonymous as they are based on 
different criteria; there exist instabilities which occur 
at large spatial scales which are driven by velocity- 
space effects. These instabilities are macroscopic by 
the Treumann definition and velocity-space instabil¬ 
ities by the Krall definition. The mirror instability is 
such an instability, which arises across a broad spec¬ 
trum of scales (see the bottom row of Figure [2]), but 
whose mechanism relies on a resonant response of a 
narrow region of the velocity distribution functioni^ 
Here we use the macroscopic/microscopic terminol¬ 
ogy solely to classify the spatial scales at which 
the instability arises and the configuration/velocity 
space terminology to classify the nature of the insta¬ 
bility mechanism. 

The double adiabatic Chew-Goldberger-Low^ 
(CGL) firehose instability is a canonical example of 
a macroscopic, configuration-space instability. The 
instability arises for large spatial scales with kpp <C 
1, as shown in row three, center column of Fig¬ 
ure m for /3|| > 1 and Tj_/T\\ < 1. The mech¬ 
anism driving the instability, that the perpendicu¬ 
lar pressure is not sufficient to counteract the cen¬ 
trifugal force from a bent magnetic flux tube, re¬ 
lies only on the macroscopic quantities of perpen¬ 
dicular and parallel pressure, and therefore depends 
on /3|| and Tl/T|| < 1 (the single-fluid analogs of 
/3||p and T±p/T\\p < 1). Though the parallel hre- 
hose and the Alfven firehose instabilities occur in 
the same region of i/3\\p,T±p/T\\p) parameter space, 
they are both fundamentally different from the CGL 
instability. Both are microscopic, velocity-space in¬ 
stabilities which depend on the resonance condition 
w(k) — fc||U|| ± Op = 0 and not on the bulk, ther¬ 
modynamic quantities of perpendicular and paral¬ 
lel pressure. The wavevectors for which the micro¬ 
scopic firehose instabilities arise are shown in the 
second and third rows of Figured Both the parallel 
and Alfven firehose instabilities have marginal sta¬ 
bility thresholds closer to temperature isotropy than 
the CGL firehose, meaning that for a plasma with 
T±p/T\\p decreasing from unity, the velocity-space 


instabilities will arise before the configuration-space 
instability. 


A. Proton Cyclotron Instability 

The proton cyclotron instability is a microscopic, 
velocity-space instability driven by the resonant con¬ 
dition a;(k) — fc||U|| — Hp = 0. This instability couples 
to the proton cyclotron wave^^*^ the extension of 
the Alfven wave at small parallel scales, fc||dp > 1. 
It occurs for Tj_p/T||p > I and for all /3||p, with the 
marginal stability threshold decreasing with increas¬ 
ing P\\p- 

The top row of Figure 0] presents the fast and 
Alfvenic roots in complex frequency space for plasma 
parameters relevant to the proton cyclotron insta¬ 
bility: /3||p = 1.0, k±pp = 10“^, fc||dp = 0.7, 
and rj_p/T||p = 1.00 (left column), 1.99 (center), 
and 2.51 (right). The slow and entropy modes are 
too heavily damped to be shown in these plots. 
For Tj_p/T||p = 1, the proton cyclotron waves are 
damped at a rate of y/Hp Ri —0.2. Once the per¬ 
pendicular temperature is equal to twice the parallel, 
the mode transitions from damped to growing, with 
increased anisotropy leading to an increased growth 
rate of y/Hp ~ 0.1 at Tj_p/T||p = 2.5. The damping 
of the fast waves is largely unaffected by the pres¬ 
ence of the anisotropy, with a slight increase in the 
real frequency of the fast waves. 

The microscopic nature of the proton cyclotron 
instability is shown in the top row of Figure [2J 
which plots the growth rate y/Hp > 0 as a function 
of wavevector for unstable modes with /3||p = 1.7 
and Tlp/Tiip = 1.344 (left column) and Tj_p/T||p 
= 2.688 (center column). These parameters, in¬ 
dicated as blue dots (right column), are compared 
to the associated instability threshold (black line) 
in the (/3||p, Tj_p/T||p) plane. The instability is con¬ 
fined to 0.1 < fc||(ip < 1 for both values of Tj_p/T||p. 
While the unstable region for the marginally un¬ 
stable case is restricted to a nearly parallel region 
with kj_pp <C 1, the unstable wavemodes for the 
highly unstable case include more oblique modes 
with k±pp < 1.0. It should be noted that the posi¬ 
tion in the (/3||p, Tlp/Tjip ) plasma parameter space 
for the latter case is not typically observed in the so¬ 
lar wind (presumably because it is highly unstable). 

The first column of Figure [S] illustrates the 
(/3||p, Tj_p/r||p) frequency and eigenfunction depen¬ 
dence of the proton cyclotron wave for k±pp = 10“^ 
and fe||(ip = 0.6. The frequency of the wave (top 
row) is fairly well constrained by O.lHp < w ^ Hp, 
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with some exceptions for portions of the j3\\p > 1, 
T^p/T\\p < 1 curves. The damping rate ^/Clp (sec¬ 
ond row) of the stable modes is largely insensitive to 
T^_p/T\\p but is strongly dependent on /3||p, while the 
transition from damping to growth is only dependent 
on rLp/T||p. The wave is left-hand circularly polar¬ 
ized with \5B^\ ~ \SBy\ and \5B^\ -C \5By\ (third 
through fifth rows). This circular polarization is ut¬ 
terly insensitive to changes in either j3\\p or rj_p/T||p 
over the wide range plotted. The mode is also nearly 
incompressible, with no significant increases in com¬ 
pressibility due to deviations from isotropy (sixth 
row). Note that the (i\\p = 10 (red), Tlp/Tjip < 1 
unstable mode (second row) is due to exceptional 
points of the type discussed in Appendix A, and 
is a manifestation of the parallel firehose instability 
rather than the proton cyclotron instability. 


B. Parallel Firehose Instability 

The parallel firehose instability is a microscopic, 
velocity-space instability driven by the resonant con¬ 
dition u;(k) — fc||U|| -|-flp = 0. This instability couples 
to the whistler wave^^i*^ the extension of the fast 
magnetosonic wave at small parallel scales, k\^dp > 1. 
The parallel firehose instability occurs for /3||p > 1 
and T^p/T||p < 1, with the marginal stability thresh¬ 
old decreasing for increasing /3||p. 

The second row of Figured] presents the fast and 
Alfven roots in complex frequency space for plasma 
parameters relevant to the parallel firehose insta¬ 
bility: Pwp = 2.8, k_\_pp = 10“^, k\\di = 0.5, and 
T^-p|T\\p = 1.00 (left column), 0.63 (center), and 
0.39 (right). As with the proton cyclotron insta¬ 
bility plots, the slow and entropy modes are too 
heavily damped to be shown. The fast mode passes 
through marginal stability at T^p/Typ = 0.63 and 
becomes more unstable for smaller Tip/Tup. The 
Alfven root represents the proton cyclotron wave, 
which becomes more strongly damped with decreas¬ 
ing Tlp/Tjip. In addition, the real frequency of 
the fast mode decreases slightly with decreasing 

TAp/Tiip. 

The microscopic nature of the parallel firehose in¬ 
stability is shown in the second row of Figure H] 
which plots the growth rate y/flp > 0 as a function 
of wavevector for unstable modes with j3\\p = 1.7 
and Tlp/Tiip = 0.55 (left column) and T^p/Typ = 
0.275 (center column). These parameters, indicated 
as blue dots (right column), are compared to the 
associated instability threshold (black line) in the 
(/I||p, Tj_p/r||p) plane. As with the proton cyclotron 


instability, unstable modes for the parallel firehose 
instability are confined to 0.1 < fey dp < 1. The per¬ 
pendicular extent of the instability does not signifi¬ 
cantly grow with smaller T^p/Typ, with a maximum 
perpendicular scale around k_\_pp ~ 0.1, yielding un¬ 
stable wavevectors within an angle 0 ~ tan“^(0.1) ~ 
6° from the equilibrium magnetic field Bg. The 
growth rate increases with decreasing Tlp/Tup, but 
it never signihcantly exceeds y/flp = 0.05. Although 
the extent of the unstable region in (fc_LPp, fc||dp) 
wavevector space for the parallel hrehose instabil¬ 
ity is similar to that of the other parallel instability, 
the proton cyclotron instability (compare first and 
second rows), for the more unstable cases (center 
column) the proton cyclotron instability extends to 
more oblique angles of 0 < 45°, rather than the limit 
of 0 < 6° for the parallel firehose instability. 

The second column of Figure [5] illustrates the 
(/S||p, Tlp/Tjip) frequency and eigenfunction depen¬ 
dence of the parallel whistler wave for k±pp = 10“^ 
and fell dp = 0.5. The frequency (top row) is very 
tightly constrained to w/Op ~ 1 regardless of /3||p or 
rLp/T||p. There is a significant dependence of y/flp 
(second row) on T^p/Typ for I3\\p > 1. The par¬ 
allel whistler waves driven by the parallel firehose 
instability are right-hand circularly polarized, with 
\5Bj:\ — \5By\ and \5Bz\ <?; \dBy\ (third through 
fifth rows), and are nearly incompressible (sixth 
row). These magnetic and density eigenfunctions of 
the whistler waves associated with the parallel fire¬ 
hose instability are insensitive to changes in j3\\p or 

T^p/T^p. 


C. Alfven Firehose Instability 

The Alfven firehose instability is a microscopic, 
velocity-space instability driven by the resonant con¬ 
dition a;(k) — A:||U|| ± flp = This instability 
couples to non-propagating, w = 0 oblique Alfven 
waves, and occurs for /3||p > 1 and T^p/T\\p < 1. 
The zero real frequency region of the Alfven disper¬ 
sion surface grows with increasing /3||p and decreas¬ 
ing Tj_p/T||p, transforming the resonant condition to 
fc|| = ±.rtp/v\\. 

The third row of Figure|4]presents the Alfven, fast, 
and entropy roots in complex frequency space for pa¬ 
rameters relevant to the Alfven firehose instability: 
,5||p = 2.8, fcj_pp = 0.45, A:||dp = 0.34, and Tlp/T||p 
= 1 (left column), 0.39 (center), and 0.25 (right). 
The real frequency of the propagating Alfven waves 
decreases for decreasing T^p/T||p until the ±a; modes 
become degenerate at a; = 0 at the marginally stable 
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state. A further decrease in T±^p/T\\p breaks the de¬ 
generacy of these two Alfven solutions, causing one 
of the modes to damp and the other to grow unsta¬ 
ble. Both stable and unstable modes remain non¬ 
propagating with real frequency w = 0. As the non¬ 
linear saturation of the instability pushes the plasma 
from an unstable point in (/?||p, rj_p/T||p) parameter 
space back towards a state of marginal stability, the 
instability-driven, non-propagating modes transition 
back into propagating oblique Alfven wavesi^ 

The wavevector regions of unstable growth for the 
Alfven and CGL firehose instabilities are shown in 
the third row of Figure [2] Plasma parameters are 
set to /3||p = 2.8 and Tlp/Typ = 0.4 (left column) 
and 'Tip/Tiip = 0.2 (center). These parameters, indi¬ 
cated as blue dots (right column), are plotted against 
the Alfven (red solid) and CGL (black dash-dot) 
firehose instability thresholds. For the marginally 
unstable plasma (left), only the microscopic Alfven 
hrehose instability arises, with the unstable region 
tightly constrained to a narrow teardrop-shaped re¬ 
gion near fcj_pp ~ k\\dp ^ 1. For the more un¬ 
stable case with T^p/T\\p = 0.2 (center), both the 
microscopic Alfven and macroscopic CGL instabil¬ 
ities occur. The CGL hrehose instability causes a 
broad spectrum of unstable modes with fc||fip < 0.1 
and k±pp < 0.6, while the Alfven hrehose is still re¬ 
stricted to the oblique wavevectors at kinetic scales 
with k_\_pp ~ A:||dp ~ 1. The CGL hrehose un¬ 
stable modes have a smaller growth rate than the 
Alfven hrehose unstable modes, but the unstable 
CGL modes are very widely distributed in wavevec¬ 
tor space, potentially allowing them to dynami¬ 
cally interact with the anisotropic huctuations of the 
large-scale turbulent cascade. 

The third column of Figure [S] illustrates the 
(/?l|p, Tj_p/T||p) frequency and eigenfunction depen¬ 
dence of the oblique kinetic scale Alfven wave for 
kj_Pp = 0.1 and k\\dp = 0.1. For j3\\p < 1, the real 
frequency (top row) uj/D,p ~ 0.1, while for increasing 
/3||p, uj/Qp increases for rj_p/T||p > 1 and decreases 
for Tj_p/T||p < 1, asymptoting to a; = 0 in the unsta¬ 
ble region. Note that for /3||p = 10 (red), the Alfven 
mode transitions between propagating (w ^ 0) and 
non-propagating (w = 0) several times as T±p/T\^p 
decreases from 1 to 0.1. The /3||p = 10, Tj_p/T||p 
> 1 unstable mode (second row) is the proton cy¬ 
clotron instability, which reaches oblique angles for 
very high degrees of anisotropy. The magnetic eigen¬ 
function (third through fifth rows) for the stable 
Alfven wave (with ,5||p < 1) is dominated by 5By, 
with a transition to elliptical polarization for the 
Alfven firehose unstable modes at /3||p > 1. While 
more compressible than the two parallel unstable 


modes, the oblique Alfven wave is still mostly incom¬ 
pressible (sixth row). The magnetic eigenfunctions 
associated with the Alfven firehose instability have 
signihcantly more dependence on /3||p and rj_p/T||p 
than those associated with the parallel instabilities. 


D. Mirror Instability 

The mirror instability is a macroscopic, velocity- 
space instability which couples to oblique entropy 
modesj ^^i^° It occurs for Tj_p/r||p > 1 and for all /3||p, 
with the marginal stability threshold decreasing for 
increasing ,5||p. 

The mirror instability was originally proposed as 
a configuration-space instability^!^ which was trig¬ 
gered by the anti-phase correlation between macro¬ 
scopic changes in the pressure and magnetic field 
strength. The instability was shown to couple with a 
non-propagating mode, though it was thought that 
the associated stable mode was propagating. This 
mode was identified as the MHD slow wave. Kinetic 
description a^^d^ have demonstrated that the physics 
of the configuration-space description was incorrect, 
due to its inability to model the velocity-space ef¬ 
fects driving the instability. Unfortunately, the con¬ 
ception of the mirror instability as a configuration- 
space instability has persisted in the literature for 
several decades. The instability threshold calculated 
from configuration-space theory is still correct, but 
the mode driven unstable is the entropy mode, not 
the slow wavei^ 

The mirror instability actually arises due to the 
difference between the anti-phase response of the 
bulk plasma’s thermal pressure to magnetic pressure 
perturbations and the in-phase response of particles 
with U|| ~ 0. The particles with ?;|| ~ 0, typically re¬ 
ferred to as resonant particles, move very little along 
the magnetic field lines, gaining or losing energy with 
increasing or decreasing field strength. Particles in 
the bulk of the plasma with finite uy stream along 
the magnetic field lines, largely conserving particle 
energy by transferring it between their parallel and 
perpendicular degrees of freedom. A through review 
and discussion of the physics relevant to the mirror 
instability can be found in Southwood & Kivelson 
1993 ,^ 

The marginal stability threshold occurs at larger 
Ti.p/T\\p for the mirror instability than for the pro¬ 
ton cyclotron instability (see Figure [T|), leading to 
the contentionii that the larger proton temperature 
anisotropies needed to drive the mirror instability 
would be isotropized by the proton cyclotron in- 
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stability before the mirror instability could be trig¬ 
gered. However, the inclusion of minor ions^ can 
alter the instability thresholds, reducing the lin¬ 
ear growth rate of the proton cyclotron instabil¬ 
ity while leaving the mirror instability mostly un¬ 
changed. Additionally, solar wind measurements of 
(/3||p, Tj_p/T||p) usually appear to be constrained by 
the marginal stability threshold of the mirror in¬ 
stability rather than that of the proton cyclotron 
instabilityj^i^ Which of these two instabilities dom¬ 
inates in the free-streaming solar wind and what 
mechanisms control their interaction remains an 
open scientific question. 

The bottom row of Figured presents the Alfven, 
fast, slow, and entropy roots in complex frequency 
space for parameters relevant to the mirror insta¬ 
bility: j3\\p = 1.0, fc_L/Op = 0.15, k\\dp = 0.4, and 
'Txp/T'iip = 1.0 (left column), 1.58 (center), and 2.51 
(right). For these parameters, the entropy mode 
is less damped than the slow modes at T^p/T\\p 
= 1.0, and as T^p/Typ increases, the damping rate 
of the entropy mode decreases to zero, reaching the 
marginally stable state at Tlp/Tjip = 1.58. Above 
this value, the mirror instability is triggered, result¬ 
ing in a growing, non-propagating (w = 0) mode. 
Like the Alfven firehose instability, the mirror in¬ 
stability passes through (a;, 7 ) = (0,0), but un¬ 
like the Alfven firehose, the related stable mode 
is non-propagating, meaning that the return to a 
marginally stable state of the plasma will not result 
in the production of propagating waves. For most 
values of T^p/Typ, the Alfven, fast, and slow modes 
are essentially unchanged, with the exception of very 
large temperature anisotropies. For T^p/Typ ~ 10, 
the proton cyclotron instability goes unstable, even 
at the oblique wavevector angles relevant to the mir¬ 
ror instability, as seen in the upper right panel of 
Figure El 

The bottom row of Figure E] illustrates the macro¬ 
scopic nature of the mirror instability by plotting the 
growth rate y/Hp > 0 as a function of wavevector for 
unstable modes with /?||p = 1.0 and T^p/Typ = 2.0 
(left column) and 4.0 (center). In both cases, the 
unstable modes fill a broad region with fey dp < k^pp 
and k±pp < 1. Within the unstable region, the 
growth rate y/Hp increases linearly with A:||dp, so 
the most rapid growth occurs for the largest unsta¬ 
ble values of fey dp. Increased anisotropy leads to an 
increase in the growth rate and increases the area of 
the unstable wavevector region slightly, though the 
unstable region remains bounded by k±pp < I. 

The fourth column of Figure [5] illustrates the 
(/?IIP ,T±p/T\\p) frequency and eigenfunction depen¬ 
dence of the entropy mode for k±pp = 0.6 


and fcydp = 0.1. The entropy modes are non¬ 
propagating, with zero real frequency a; = 0, for 
all /3||p and Tj_p/T[|p (top row). The damping rate 
(second row) has a slight dependence on /3||p, with 
larger /3||p generally leading to larger damping or 
growth rates. For a fixed wavevector, the anisotropy 
value rLp/Tjip at marginal stability increases with 
decreasing /3||p. The magnetic eigenfunctions for the 
mirror-instability-driven entropy modes are gener¬ 
ally dominated by with a subdominant con¬ 

tribution from SBj; which increases with increasing 
/3||p. The mirror-instability-driven entropy mode is 
the most compressible of the four modes considered 
in this investigation, with larger density fluctuations 
for smaller /3||p and larger Tj_p/T||p. 

IV. IMPACT OF PROTON 
TEMPERATURE ANISOTROPY ON THE 
LARGE-SCALE CASCADE 

Now that we have established the properties of 
the four proton temperature anisotropy instabilities, 
we consider how the proton temperature anisotropy 
impacts the nonlinear dynamics of the Alfvenic fluc¬ 
tuations underlying the large-scale cascade. 

Early research on incompressible MHD turbulence 
in the I960 a^^'^^ emphasized the wave-like nature of 
turbulent plasma motions, suggesting that nonlin¬ 
ear interactions between counterpropagating Alfven 
waves—or Alfven wave collisions—mediate the tur¬ 
bulent cascade of energy from large to small scales. 
Following significant previous studies on weak in¬ 
compressible MHD turbulence the nonlinear 

energy transfer in Alfven wave collisions has re¬ 
cently been solved analytically in the weakly non¬ 
linear limit confirmed numerically with gyroki- 
netic simulations?^ and verified experimentally in 
the laboratory?^^— establishing Alfven wave colli¬ 
sions as the fundamental building block of astro- 
physical plasma turbulence. A discussion of the im¬ 
portant role of linear wave physics in the context of 
strong turbulence can be found in Howes, Klein, & 
TenBarge 2014^2 and Howes 2015.-^ 

The physics of Alfven wave collisions in the iner¬ 
tial range is most clearly illustrated in the context 
of the incompressible MHD equations expressed in 
the symmetrized Elsasser form?^ 

Bz^ 

— T • Vz± = -z^ • Vz± - VP/po, (4) 

and V • z^ = 0. Here = Bo/-\/47rpo is the Alfven 
velocity due to the equilibrium field Bg = BqZ where 
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B = Bo + (5B, P is total pressure (thermal plus mag¬ 
netic), po is mass density, and = u± S’B/^yAnpo 
are the Elsasser fields which represent waves that 
propagate up or down the mean magnetic field. The 
• Vz^ term governs the nonlinear interactions be¬ 
tween counterpropagating Alfven waves, or Alfven 
wave collisions. 

The mathematical form of the nonlinear term, 
z^ • Vz^, determines important properties that gov¬ 
ern the turbulent cascade of energy from large to 
small scales. The first, and most fundamental, prop¬ 
erty is that only counterpropagating waves inter¬ 
act nonlinearly )^^i^^i^^i^° a property that is possi¬ 
ble because Alfven waves (and pseudo-Alfven waves 
in the case of incompressible MHD) are not disper¬ 
sive. If not, waves traveling in the same direction 
could catch up with each other and interact nonlin- 
early. Second, the property that only counterprop¬ 
agating waves interact nonlinearly fundamentally 
leads to an anisotropic cascade of energy in plasma 
turbulence?^ in which energy is preferentially trans¬ 
ferred to small perpendicular scales, leading to 
small-scale turbulent fluctuations with an anisotropy 
k± ^ fc||. In this anisotropic limit, the Alfven waves 
dominate the nonlinear interactions?^!^^— with 
compressible fluctuations relegated to a subdomi¬ 
nant role in the nonlinear energy transfer. Third, 
in order for the nonlinearity to be nonzero, the vec¬ 
tor nature of the nonlinear term requires that two 
counterpropagating Alfven waves, each of which is 
linearly polarized in the MHD limit kpp -C I, not be 
polarized in the same planei^ Therefore, to under¬ 
stand how the proton temperature anisotropy im¬ 
pacts the large-scale turbulent cascade, it is impor¬ 
tant to explore how Tj_p/T||p 1 affects these two 
important properties of the anisotropic Alfven waves 
that control the turbulent cascade: (I) the nondis- 
persive nature of the Alfven waves, and (2) the po¬ 
larization of the Alfven wave eigenfunction. 

The large-scale cascade is dominated by Alfvenic 
fluctuations in an anisotropic region of wavevector 
space with k± ^ fc||. This sense of anisotropy 
is well-supported by multi-spacecraft measurements 
of turbulence near the ion kinetic scales in the so¬ 
lar windj^^— Modern scaling theories of Alfvenic 
turbulence suggest that a critical balance between 
linear and nonlinear timescales leads to turbulent 
fluctuations that obey a scale-dependent wavevec¬ 
tor anisotropy . To explore the effect of the 
proton temperature anisotropy Tj_p/T||p on these 
anisotropic Alfven waves, we solve for their proper¬ 
ties along an idealized line of critical balance given 


byl^~— 


^\\Pp ~ i^OPp) ^ 


(k\_Ppfl^ + {k±ppy/^ ' 
I + {k^PpY 


(5) 


where k^pp = 10“^ is the isotropic driving scale 
(the outer scale of the turbulent inertial range), at 
which fluctuations are assumed to be isotropic with 
fcj- = ^11 ■ This critical balance relation represents 
the approximate upper boundary (in fc||) in wavevec¬ 
tor space of the power in the anisotropic turbu¬ 
lent fluctuations , Limits of this idealized relation 

lead to a wavevector anisotropy scaling as k\\ oc 
in the MHD inertial range {k±pp <C 1) and A:|| oc 
in the kinetic dissipation range {k±pp > 1)4^ The 
wavevector scalings in those two limits are supported 
by MHD simulations of Alfven wave turbulenc o^^i^^ 
and gyrokinetic simulations of kinetic Alfven wave 
turbulence!^ 

In Figure I we plot the real frequency as well 
as two metrics of the eigenfunction polarization for 
Alfven waves for the Vlasov-Maxwell system along 
this idealized critical-balance line. The columns are 
organized by /3||p, with /3||p = 0.1 (left column), /3||p 
= 1.0 (center), and /3||p = 10.0 (right). We explore 
the range of proton temperature anisotropies that 
are observed in the solar wind^ at each value of /3||p, 
so each column has different values for {T^p/T\\p)yj^\y^ 
and (T'j_p/r||p)max. The temperature anisotropy is 
indicated by the color of each curve, with the color- 
bar displaying the variation of colors from the mini¬ 
mum to the maximum temperature anisotropy for 
that column. The minimum proton temperature 
anisotropy (r'-Lp/T||p)jjjjjj is blue, an isotropic pro¬ 
ton temperature Tj_p/T||p = 1 is black, and the max¬ 
imum (T^p/T'||p)max is red. 

Alfven waves propagate along the local mean mag¬ 
netic field Bq at the parallel group velocity Vg\\ = 
duj/dk^^. If this group velocity is independent of the 
wavevector, then the Alfven waves are nondispersive 
and an arbitrary Alfven wavepacket will propagate 
without distortion along the magnetic field at the 
parallel group velocity. In the top row of Figure [6l 
the normalized frequency uj/k^^VA is plotted vs. nor¬ 
malized perpendicular wavenumber kj_pp along the 
path of critical balance given by ®. For each value 
of /?||p, we see that the Alfven waves are indeed 
nondispersive (the normalized curve uj/k^^VA is inde¬ 
pendent of kj_, and therefore Ug|| is constant), except 
for k±pp > 1 where finite proton Larmor radius ef¬ 
fects lead to a transition to dispersive kinetic Alfven 
waves. 

For an isotropic proton temperature T±p/T\\p = I, 
the Alfven waves propagate at the Alfven velocity 
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VA for any value of (3\\p (black). For an anisotropic 
proton temperature T^p/Typ ^ 1, the parallel group 
velocity depends on the value of T±^p/T\\p (partic¬ 
ularly for ^iip > 1), but for any particular choice 
of plasma parameters (/3||p,Tj_p/T||p), the Alfven 
wave remains nondispersive in the inertial range 
fc±Pp ^ 1- In fact, the dispersion relation is well 
modeled by the double adiabatic CGL resultj^i^ 
w = k\\VA^J^ + /3||(7j_/T|| — l)/2. Therefore, the 
nondispersive nature of Alfven waves is insensitive 
to the proton temperature anisotropy, and so we ex¬ 
pect that the property that only counterpropagating 
waves interact nonlinearly persists for Tip/Tjip ^ 1. 
Of course, when Alfven waves become unstable to 
the macroscopic CGL firehose instability, they be¬ 
come non-propagating with w = 0, as can be seen 
for the two lowest values of Tj_p/T||p (blue) in the 
f5\\p = 10 column of Figure [H However, a nearly 
negligible fraction of observed solar wind intervals 
have parameters unstable to the CGL firehose^^— 
reducing the significance of these non-propagating 
modes for the large-scale turbulent cascade. It is 
worthwhile emphasizing, also, that the parallel fire¬ 
hose and Alfven firehose velocity-space instabilities 
are unstable in regions of wavevector space (see Fig¬ 
ure El) that do not overlap the anisotropic region of 
wavevector space inhabited by the Alfvenic fluctua¬ 
tions associated with the large-scale cascade, given 
by modes with parallel wavenumbers below the crit¬ 
ical balance line ([5]). 

Next, we investigate whether an anisotropic pro¬ 
ton temperature alters the polarization of the 
anisotropic Alfven waves associated with the large- 
scale cascade. For an isotropic temperature T^p/Tup 
= 1, Alfven waves are linearly polarized, with 
\6By\ ^ \^Bx\ ~ \^Bz\. To determine any 

anisotropy-induced deviations from linear polariza¬ 
tion, we plot \5Bx\/\5By\ in the second row of 
Figure [6] for the same range of plasma parame¬ 
ters (/3||p,Tj_p/r||p). The Alfven wave remains lin¬ 
early polarized, with \5By\ ^ \5Bx\^ with little 
change from the isotropic temperature case T^p/T\\p 
= I (black) for all values in the observed range of 
(/3||p,T'j-p/r'||p)- This result holds even for the modes 
unstable to the GGL firehose in the /3||p = 10 case. 
With no changes in the polarization of the Alfven 
waves, we expect no significant modifications of the 
nonlinearity responsible for the large-scale turbulent 
cascade due to proton temperature anisotropy. 

Finally, for an isotropic temperature Tj_p/T||p = 1, 
the eigenfunction relation for Alfven waves is given 
by SBy/Bo = ±SuylvA, where the sign dictates the 
direction of propagation of the Alfven wave along 
the mean magnetic field Bg. This property has been 


used to identify large-scale Alfven waves in the solar 
wind^ and is the physical basis enabling upward and 
downward propagating Alfven waves to be described 
by the Elsasser variables, jvA = vl/va ± 

Using the double adiabatic GGL dispersion rela¬ 
tion for MHD Alfven waves, we obtain an eigen¬ 
value relation generalized to account for tempera¬ 
ture anisotropy, to/{k\\VA)5By/B q = ±5uylvA- This 
simple result enables one to construct generalized El¬ 
sasser variables for the case of anisotropic tem¬ 
perature, ZqIva = u/u^ ± a;/(/c||UA)(5B/Hg. We 
test this generalization of the Elsasser variables by 
plotting a;/(fc||ri^)|(5i?p/Hg|/|i5np/uA| vs. fc^pp in the 
third row of Figure |6l demonstrating clearly the 
validity of this generalization of the Elsasser vari¬ 
ables for Alfven waves throughout the inertial range 
fc_LPp ^ 1 for any choice of plasma parameters 

(/3||p,Tlp/T||p). 


In summary, we find that neither the nondisper¬ 
sive nature nor the polarization of the Alfven waves 
that constitute the large-scale turbulent cascade 
are altered by the proton temperature anisotropy. 
Therefore, we conclude that the physics of the large- 
scale cascade in the inertial range is insensitive to 
the proton temperature anisotropy over the range of 
values observed in the solar wind. Gonsequently, we 
expect that studies of Alfven wave nonlinear turbu¬ 
lent interactions using an isotropic proton velocity 
distribution will still accurately describe turbulence 
in the inertial range properly, even for plasmas with 

Tlp/T||p ^ 1. 


The empirical prediction presented in this section, 
that the turbulent dynamics of the large-scale cas¬ 
cade is not significantly altered by proton temper¬ 
ature anisotropy, is supported by the results of a 
recent theoretical treatment of kinetic turbulence 
in the inertial range for plasmas with anisotropic 
temperature distributions or relative drift among 
ion species^i In that study, the main physical fea¬ 
tures of plasma turbulence in the inertial range per¬ 
sist for non-Maxwellian distributions: the Alfvenic 
and compressive fluctuations are decoupled, with 
the latter passively advected by the former; and the 
Alfvenic cascade remains essentially fluid with dy¬ 
namics governed by the equations of reduced MHD 
modified to account for changes in the Alfven speed 
due to pressure anisotropy and species drifts. 
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V. ENERGY INJECTION BY PROTON 
TEMPERATURE ANISOTROPY 
INSTABILITIES 

Although the proton temperature anisotropy does 
not appear to alter significantly the nature of the 
turbulent nonlinear interactions involved in the 
anisotropic cascade of energy from large to small 
scales, kinetic instabilities can drive electromagnetic 
fluctuations in the plasma, thereby directly inject¬ 
ing energy into fluctuations at kinetic scales. Here 
we address two key questions about these instability- 
driven fluctuations: (1) What are the observable sig¬ 
natures of the instability-driven fluctuations, espe¬ 
cially considering that spacecraft measurements are 
made both in a moving frame of reference and in 
the presence of the large-scale turbulent cascade?; 
and (2) How do these instability-driven fluctuations 
interact nonlinearly with each other and with the 
fluctuations of the large-scale cascade? 


A. Observable Signatures 

Let us consider first the instability-driven fluctua¬ 
tions in the frame of the solar wind plasma, a frame 
of reference generally moving at a super-Alfvenic ve¬ 
locity with respect to the spacecraft frame in which 
measurements are madei ^^’^^ The two parallel pro¬ 
ton temperature anisotropy instabilities, the proton 
cyclotron and parallel firehose instabilities, gener¬ 
ate propagating proton cyclotron and whistler waves 
with nonzero frequencies (see lower panels in Fig¬ 
ure [a in the wavevector region with 0.1 < /c||dp < 1 
and k± <g; fc|| (see top two rows of Figure [2]). These 
fluctuations may serve to transport thermal energy 
from one spatial region of the solar wind to an¬ 
other via the Poynting flux of the unstable waves. 
On the other hand, the two oblique instabilities, 
the Alfven firehose and mirror instabilities, gener¬ 
ate non-propagating fluctuations with a; = 0 in the 
plasma frame (see upper panels of Figure [S]). The 
Alfven firehose instability drives fluctuations with 
peak growth rates in the wavevector region with 
0.1 < fc||(ip < 0.4 and 0.1 < kj_pp < 1 (see third row 
of Figure [2]). The mirror instability drives fluctua¬ 
tions over a broad wavevector region with fey < k± 
and k±pp ;< 1, but its peak growth rates (which in¬ 
crease linearly with A:||) are restricted to the small¬ 
est scales, 0.1 < k±pp < 1 (see bottom row of Fig¬ 
ure [2]). We expect that the more slowly growing mir¬ 
ror modes at k±pp < 0.1 will not be observable be¬ 
cause the fluctuations of the large-scale cascade oc¬ 


cupying the same region of wavevector space, which 
increase in amplitude with decreasing wavenumber, 
will be significantly larger and therefore dominate 
measurements at these low perpendicular wavenum¬ 
bers. 

Next we consider the signature of these instability- 
driven fluctuations as measured in the spacecraft 
frame as the solar wind plasma flows past at a super- 
Alfvenic velocity. In general, for a fluctuation with 
wavevector k and plasma-frame frequency oj flow¬ 
ing past the sampling spacecraft at solar wind ve¬ 
locity vsw, the spacecraft-frame frequency is given 
by (jJsc = w -f k • Vsm, the sum of the plasma-frame 
frequency plus a Doppler-shifted spatial variation4S 
For a sufhciently fast flow speed, the plasma-frame 
frequency term is negligible compared to the Doppler 
shift term, so lOsc — k-v^u,, an approximation known 
as the Taylor hypothesis)^ widely employed in the 
solar wind since the solar wind speed is typically 
Vsw ~ lOwA- The condition for the Taylor hy¬ 
pothesis to be valid for both Alfven and whistler 
waves^ is {vsw/va) cosOkv ^ k^^dp, where k • Vgw = 
kvsw cos 9kv Since all of these instability-driven 
modes satisfy k\\dp ^ 1, we may safely adopt the 
Taylor hypothesis. An important point to empha¬ 
size here is that for choices of angle 9kv such that 
cos 9kv 1 (a condition that would violate the Tay¬ 
lor hypothesis), the spacecraft-frame frequency of 
the instability-driven fluctuation is downshifted suf¬ 
ficiently that it will be unmeasurable in the presence 
of the large-scale cascade. 

If we express the wavevector of a particular 
instability-driven mode as k = fc||b -f k^, where the 
direction of the local mean magnetic field is given 
by the unit vector b = Bq/Hq, the spacecraft-frame 
frequency can be written 

ujsc — k\\ Vsw COS 9yB + • ^rsw , (6) 

where the angle is given by Vsw ■ Bq = cos 9yB ■ 
For the proton cyclotron and parallel firehose in¬ 
stabilities, the wavevectors of unstable modes sat¬ 
isfy k± A:||, so this expression reduces to ujsc — 
k^^Vsw cos 9yB- Note that the solar wind speed Vsw 
and the angle between the solar wind flow and the 
local mean magnetic field 9yB are both directly mea¬ 
surable from spacecraft measurements, so this cal¬ 
culation establishes a direct relation between the 
spacecraft-frame frequency and parallel wavenumber 
for the parallel modes arising from these two insta¬ 
bilities. 

We must also consider whether these fluctuations 
can be measured in the presence of the broadband 
fluctuations due to the large-scale cascade. The 
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amplitude of the background turbulent fluctuations 
decreases with increasing frequency, so these par¬ 
allel instability-driven waves are most likely to be 
measurable at their maximum spacecraft-frame fre¬ 
quency ujsc — k^^sw, occurring when O^b —t 0, or 
physically when the direction of the local mean mag¬ 
netic field aligns with the solar wind flow. Such 
alignment between the magnetic field and the solar 
wind flow, OyB 0, has previously been exploited to 
study the variation of solar wind turbulence parallel 
to the magnetic field in a number of observational 
studies^^TiS^ We can use this maximum-frequency 
relation and the wavevector region of unstable modes 
to predict the linear spacecraft-frame frequency for 
the fluctuations driven by the proton cyclotron and 
parallel firehose instabilities, 

0.02^ < fsc < 0.2^, (7) 

(2p dp 

where uJsc = 

Although the Alfven firehose and mirror instabil¬ 
ities generate unstable non-propagating fluctuations 
with real frequency a; = 0 in the plasma frame, in 
the spacecraft frame these fluctuations have ujsc ^ 0 
due to the Doppler shift term. Thus, it is difficult, 
from single-point spacecraft measurements alone, to 
distinguish these non-propagating instability-driven 
fluctuations from propagating waves in the solar 
wind. The Taylor hypothesis is, of course, triv¬ 
ially satisfied for these unstable modes since a; = 0, 
yielding the equality Wsc = k • Vg^,. Again, in the 
presence of the large-scale cascade, these instability- 
driven fluctuations are most easily measured at the 
maximum spacecraft-frame frequency — kvgw 
occurring when —>■ 0. Although 9kv is not 

an observationally accessible quantity using single¬ 
point spacecraft measurements, this simple expres¬ 
sion provides a valuable estimate for the predicted 
spacecraft-frame frequency of potentially observable 
fluctuations driven by the Alfven firehose and mirror 
instabilities. Since these modes have peak growth 
rates with 0.3fe_LPp < k\\dp < k±pp for /3||p ^ 1 (see 
bottom two rows of Figure [2]), the magnitude of the 
wavevector can be estimated by fc = (fcj| -I- 
kj_, yielding a simplified prediction for the linear 
spacecraft-frame frequency of fluctuations driven by 
the Alfven firehose and mirror instabilities, 

0 . 02 ^ < < 0 . 2 ^. ( 8 ) 

Pp Pp 

Note that this expression is sensitive to the min¬ 
imum value of 9kv at which there is energy from 
instability-driven fluctuations; for example, if the so¬ 
lar wind flow is instantaneously along the local mean 


magnetic field = 0), then the unstable regions 
of wavevector space in the bottom row of Figure [5] 
would yield 9kv > Tr/d, reducing the measured uJsc 
by a factor cos9kv =0.71 from this simple estimate. 

From this analysis, a key qualitative difference be¬ 
tween the turbulent fluctuations of the large-scale 
cascade and the instability-driven fluctuations with 
solar wind relevant plasma parameters is their extent 
in frequency. The cascade of turbulent fluctuation 
energy from large to small scales is characterized by 
a broadband frequency spectrum of turbulent fluc¬ 
tuations. In contrast, the instability-driven fluctu¬ 
ations are likely to exhibit a narrowband frequency 
spectrum within the ranges given by 0 and (HI). 

Of course, whether the fluctuations driven by any 
of the four proton temperature anisotropy insta¬ 
bilities is observable in the presence of the large- 
scale cascade depends on the amplitudes of both 
the instability-driven fluctuations and the turbulent 
fluctuations comprising the large-scale cascade. The 
amplitude of the instability-driven fluctuations de¬ 
pends on the nonlinear saturation mechanism for 
each instability as well as the large-scale dynamics 
pushing the plasma into an unstable region of plasma 
parameter space, both topics beyond the scope of 
this study. 

Even without a prediction for the saturated am¬ 
plitudes of the instability-driven modes, a definitive 
identification of these modes may be possible us¬ 
ing the properties of these modes that observation- 
ally distinguish them from the dominant Alfvenic 
fluctuations and subdominant compressible fluctu¬ 
ations of the large-scale cascade. A careful com¬ 
parison of the correlations between two components 
of the measurements with that of the linear eigen¬ 
functions for each instability-driven mode can be 
used to make such an identificationiS There is a 
well-established precedent for such identifications 
using various polarizations and correlation s^^ 
in the magnetosphere22ii££ and in the bulk of the so¬ 
lar windjii^i2£ii2id£^ An example of discriminating 
properties from the eigenfunctions of the instability- 
driven waves in Figure [5] include the fact that the 
proton cyclotron waves driven by the proton cy¬ 
clotron instability have left-handed circular mag¬ 
netic polarization and are largely incompressible 
(first column), while the whistler waves driven by the 
parallel firehose instability have right-handed circu¬ 
lar magnetic polarization and are also largely in¬ 
compressible (second column). It is worthwhile not¬ 
ing that the properties of both the proton cyclotron 
waves and the whistler waves driven by these two in¬ 
stabilities are largely insensitive to both the parallel 
proton plasma beta /3||p and the proton tempera- 
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ture anisotropy Tlp/Tup, potentially making eigen¬ 
function tests of these modes particularly robust to 
plasma parameter variations. The properties of the 
fluctuations arising from the Alfven firehose and mir¬ 
ror instabilities have substantially more variation 
with changes in the plasma parameters, but it may 
be possible to devise discriminating tests for these 
modes as well. 

We conclude this section with a brief review of 
potential observations of fluctuations caused by pro¬ 
ton temperature anisotropy instabilities in the liter¬ 
ature. A wavelet analysis of measurements of the 
solar wind magnetic energy spectrum as a function 
of 9yB within high-speed streams in the ecliptic so¬ 
lar wind using the Stereo spacecraft, presented in 
Podesta 2009 shows a bump in the nearly parallel 
bin 9yB = 3° at fsc — 0.4 Hz that may be a signature 
of fluctuations driven by one of these instabilities, 
although the paper does not summarize the plasma 
parameters of the measured intervals to test whether 
they agree with the quantitative predictions of fsc 
presented here. The same wavelet analysis technique 
applied to Ulysses measurements from the fast polar 
solar wind in Wicks et al. 201022 found a bump in 
the parallel spectrum, binned over 0° < 9yB < 10°, 
at a value of kpp ~ 0.6 that may also be due to 
instability-driven fluctuations. Measurements of the 
reduced fluctuating magnetic helicityi2iii22ii£2r2i2 
sorted by period T and angle 9yB find a broad per¬ 
pendicular signature of positive magnetic helicity 
and a narrow parallel signature of negative mag¬ 
netic helicityi2iii22 The broad perpendicular signa¬ 
ture was interpreted to be due to an anisotropic dis¬ 
tribution of KAWs with fcj_ ^ fc|| associated with 
the large-scale cascade while the parallel sig¬ 
nature was proposed to arise from proton cyclotron 
waves propagating anti-sunward along the magnetic 
field2ii22 driven by the proton cyclotron instabil¬ 
ity or from whistler waves propagating sunwarcl22i22 
driven by the parallel firehose instability. A compar¬ 
ison of the magnetic helicity from these spacecraft 
measurements to that produced using the synthetic 
spacecraft data method^ has constrained the energy 
content of the instability-driven parallel modes to 
be around 5% of the energy in the large-scale cas¬ 
cade over the narrow frequency band of the paral¬ 
lel magnetic helicity signaturej222 Another example, 
although driven by a kinetic electron rather than 
a kinetic proton instability, is the recent study by 
Lacombe et al. 2014^22 which presents evidence of 
power in intermittent whistler waves at frequencies 
10 Hz < fsc ^ 100 Hz, interpreted to be driven 
by the whistler heat flux instabilityi222 These stud¬ 


ies provide significant motivation for more thorough 
investigations to identify fluctuations driven by pro¬ 
ton temperature anisotropy instabilities in the solar 
wind and to constrain the energy content of these 
fluctuations and their effect on the thermodynamic 
evolution of the solar wind. 

An important open question is how the presence 
of pre-existing turbulence impacts the linear fre¬ 
quencies and growth rates of kinetic temperature 
anisotropy instabilities. The treatment presented 
here assumes that, to lowest order, linear theory pre¬ 
dictions for the growth rates in a quiescent medium 
give reasonable values for estimation. The observa¬ 
tional evidence for modes driven by parallel insta¬ 
bilities, reviewed above, suggests that these insta¬ 
bilities do indeed operate within the turbulent solar 
wind plasma. 


B. Nonlinear Interactions 

Next we consider how the fluctuations driven by 
proton temperature anisotropy instabilities interact 
nonlinearly with each other and with the anisotropic 
fluctuations of the large-scale turbulent cascade. 

There are, in fact, several lines of reasoning that 
suggest that the instability-driven fluctuations ex¬ 
perience negligible nonlinear interactions with each 
other. First, the parallel signature seen in mag¬ 
netic helicity measurements2ii22 is found to have 
a normalized value near 1.0, suggesting that the 
instability-driven waves, if they are responsible for 
the signature, propagate almost entirely in one 
direction along the magnetic field, either anti- 
sunward for proton cyclotron waves or sunward 
for whistler waves i22i2ii22ii22 For the non-dispersive 
Alfven waves dominating the inertial range of so¬ 
lar wind turbulence, the nonlinear interactions that 
mediate the turbulent cascade of energy occur only 
between counterpropagating wavesi22“— Although 
this property does not strictly hold true for dis¬ 
persive waves (since waves with higher perpendic¬ 
ular wavenumbers can overtake waves with lower 
perpendicular wavenumbers), any nonlinear interac¬ 
tions that do occur may be weaker by virtue of the 
longer timescales for waves to overtake each other. 
Second, the energy content of the instability-driven 
waves has been found in a recent study to be only 5% 
of the large-scale cascade fluctuation energy over the 
same spacecraft-frame frequency band .222 There¬ 
fore, these unstable waves either have a small ampli¬ 
tude or are spatially intermittent. Small amplitude 
implies a weaker nonlinearity, and spatial intermit- 
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tency implies lower probability that two wavepackets 
will collide and interact nonlinearly; either way, the 
nonlinear interactions are likely to be weaker than 
those associated with the large-scale cascade. Fi¬ 
nally, the potential observations of instability-driven 
waves reviewed above are all qualitatively narrow- 
band in spacecraft-frame frequency, whereas a typ¬ 
ical characteristic of turbulence with strong nonlin¬ 
earities is a broadband spectrum of fluctuations, sug¬ 
gesting that these instability-driven waves do not un¬ 
dergo a significant turbulent cascade. 

Finally, we consider nonlinear interactions of the 
instability-driven fluctuations with the turbulent 
fluctuations of the large-scale cascade. The first key 
point is that, with the exception of the mirror in¬ 
stability, the unstable fluctuations generated by the 
proton temperature anisotropy instabilities arise in 
a region of wavevector space (generally 0.1 < k\\dp < 
1 and fc_L ^ fc||) that does not overlap with the 
anisotropic fluctuations of the large-scale cascade in 
the region k± ^ fc||. Even for the mirror instability, 
the peak growth rates occur at moderately oblique 
wavevectors with k^pp ~ 3fc||(ip, wavevectors adja¬ 
cent to, but not generally within, the region k±^ ^ k\\ 
of the large-scale cascade. 

For the propagating waves driven by the pro¬ 
ton cyclotron and parallel firehose instabilities, this 
mismatch in wavevector space leads to signihcantly 
higher plasma-frame frequencies for the instability- 
driven waves than the anisotropic, low-frequency 
fluctuations of the large-scale cascade. In Fig¬ 
ure [71 we plot the normalized plasma-frame fre¬ 
quency uj/VLp for both Alfven waves along criti¬ 
cal balance (dash-dotted lines) and the unstable 
modes generated by the parallel instabilities with 
q/rip > 10“^ (shaded region). The proton tem¬ 
perature anisotropy is varied from Tj_p/T||p = 2.68 
(red) to Tj_p/T||p = 1.00 (black) to T^p/T\\p = 0.27 
(blue), with (dy = 1.7. This significant difference 
in plasma-frame frequencies results in an effective 
impedance mismatch that is expected to prevent 
significant nonlinear interactions between the high- 
frequency unstable waves and the low-frequency tur¬ 
bulent fluctuations. Physically, high-frequency fluc¬ 
tuations experience much lower frequency fluctua¬ 
tions as a nearly static background variation, while 
low-frequency fluctuations are only weakly affected 
by the rapid oscillation of much higher frequency 
fluctuations. But, as these instability-driven waves 
attempt to propagate along the tangled magnetic 
held associated with the large-scale turbulent cas¬ 
cade, wavepackets may be sheared out, effectively 
transferring their energy to smaller perpendicular 
scales. Thus, the energy of the instability-driven 


waves may be cascaded to smaller scales by the fluc¬ 
tuations of the large-scale cascade, but the unstable 
waves are not likely to alter signihcantly the turbu¬ 
lent nonlinear dynamics of the large-scale cascade, 
especially if the instability-driven waves have small 
amplitudes or are spatially intermittent. 

Of course, the Alfven hrehose and mirror insta¬ 
bilities generate non-propagating huctuations with 
zero plasma-frame frequency, a; = 0. These huc¬ 
tuations will essentially form a static background 
upon which the active large-scale turbulent cascade 
must proceed. It is possible that these modes can al¬ 
ter the nonlinear energy transfer of the anisotropic 
Alfvenic huctuations, meriting further considera¬ 
tion, although if the growth of these unstable huc¬ 
tuations ceases at small amplitude, their ehect on 
the large-scale cascade may be minimal. On the 
other hand, these static background modes may in¬ 
deed be cascaded to smaller scales by the active 
Alfvenic huctuations of the large-scale cascade. Fu¬ 
ture study into the ehects of the turbulent cascade 
on the nonlinear saturation of these instabilities will 
be necessary to determine if the associated unstable 
modes and nonlinear structures can grow to suffi¬ 
ciently large amplitudes to ahect the large-scale cas¬ 
cade. 


VI. CONCLUSION 

In the weakly collisional solar wind plasma, non- 
Maxwellian particle velocity distributions are rou¬ 
tinely measured, but how these non-Maxwellian dis¬ 
tributions impact the physics of plasma turbulence 
in the solar wind remains unanswered. Here we take 
a hrst step in addressing this question by exploring 
in detail the ehect of a bi-Maxwellian proton temper¬ 
ature distribution on the plasma turbulence. Specih- 
cally, we aim to understand how the proton temper¬ 
ature anisotropy ahects the nonlinear dynamics of 
the large-scale turbulent cascade and how the huctu¬ 
ations driven by the proton temperature anisotropy 
interact nonlinearly with each other and with the 
huctuations of the large-scale cascade. 

For a plasma with a bi-Maxwellian proton distri¬ 
bution, there exist at least four proton temperature 
anisotropy instabilities: the proton cyclotron, paral¬ 
lel hrehose, Alfven hrehose, and mirror instabilities. 
We present a unihed framework for these instabil¬ 
ities, identifying the associated stable eigenmodes, 
highlighting the unstable region of wavevector space, 
and presenting the properties of the growing eigen¬ 
functions. 
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We find that the proton temperature anisotropy 
significantly affects neither the nondispersive na¬ 
ture nor the polarization of the Alfven waves that 
constitute the large-scale turbulent cascade, so we 
conclude that the nonlinear interactions governing 
the turbulent cascade of energy through the iner¬ 
tial range are insensitive to the proton temperature 
anisotropy. Consequently, we expect that studies of 
the nonlinear dynamics of plasma turbulence assum¬ 
ing an isotropic proton velocity distribution will still 
qualitatively and quantitatively describe the lowest- 
order behavior of solar wind turbulence in the iner¬ 
tial range, even for plasmas with T±^p/T\\p ^ 1. The 
electromagnetic fluctuations driven by these tem¬ 
perature anisotropy instabilities are expected to be 
observed within a narrowband of linear spacecraft- 
frame frequencies Q.02vsw/'rp < fsc < ^■‘^-Vsw/fp, 
where Vp is either the proton inertial length or pro¬ 
ton Larmor radius, and the properties of the unsta¬ 
ble linear eigenfunctions may be exploited to iden¬ 
tify these instability-driven fluctuations in space¬ 
craft measurements. We argue that the instability- 
driven waves are unlikely to interact nonlinearly 
with each other, and that, although their energy 
may be cascaded to smaller scale by the turbulent 
fluctuations of the large-scale cascade, it is unlikely 
that the instability-driven waves significantly alter 
the nonlinear dynamics of the large-scale cascade. 

These findings are fully consistent with a recent 
perspective on the fluid and kinetic aspects of the 
weakly collisional plasma turbulence in the solar 
wind^ as well as an analytic treatment of tempera¬ 
ture anisotropic inertial range kinetic turbulence^^ 
In this picture, the nonlinear interactions respon¬ 
sible for the turbulent cascade of energy and the 
formation of current sheets in kinetic plasma tur¬ 
bulence are essentially fluid in nature, and are 
not significantly impacted by the typically non- 
Maxwellian form of the distribution functions. On 
the other hand, the damping of the turbulent fluc¬ 
tuations of the large-scale cascade and the injec¬ 
tion of energy by kinetic instabilities via collision¬ 
less wave-particle interactions are essentially kinetic 
in nature, and thus may depend sensitively on the 
non-Maxwellian form of the distribution functions. 
This viewpoint strongly contradicts another recent 
study using a hybrid Vlasov-Maxwell numerical ap¬ 
proach which claims that the non-Maxwellian fea¬ 
tures of the proton distribution function “may be a 
key point for understanding the complex nature of 
plasma turbulence”, and that approaches employ¬ 
ing Maxwellian or bi-Maxwellian distribution func¬ 
tions have limited applicability for the study of solar 
wind turbulencefiii The study of kinetic turbulence, 


a new frontier in heliophysics research, will clearly 
remain a hot topic for the foreseeable future. 
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Appendix A: Exceptional Points and Mode 
Identification 

The delineation between distinct dispersion sur¬ 
faces of the Vlasov-Maxwell system breaks down in 
parameter regimes where the dispersion relation con¬ 
tain branch point singularities, known in this context 
as exceptional points4^ Following a path in param¬ 
eter space which encircles such an exceptional point 
results in a continuous transition between distinct 
solutions. For instance, consider following a partic¬ 
ular solution Wo along a path connecting 
{A 2 ,Bi), (^ 2 ,^ 2 ), {Ai,B 2 ), and (Ai,Ri), where A 
and B are parameters of the Vlasov-Maxwell sys¬ 
tem and Ai _2 and B 12 are chosen to enclose the 
exceptional point {A*, ’b*); A* G (^ 1 ,^ 2 ) and B* € 
(Bi,B 2 ). The solution wi found upon returning to 
(Ai, Bi) will not be the same as the initial solution 
wq. To return to wq requires multiple revolutions 
around the exceptional point. In general, excep¬ 
tional points signify that the dispersion surfaces are 
not disjoint, but that the surfaces are multi-valued 
functions of the system parameters. 

Exceptional points for a particular system occur 
when the roots of the system, defined by the equa¬ 
tion \'D{uj)\ = 0 also satisfy ^|T>(w)| = 0,— where V 
is an n X n matrix defined by the system’s equations 
of motion. In the case of a Vlasov-Maxwell system 
|2?(w)| is the dispersion relation derived from the 
homogeneous wave equation, given by Equation 10- 
73 in Stix 1992.— As the Vlasov-Maxwell system of 
equations is complicated, we leave numerical or an¬ 
alytical treatments for the locations of exceptional 
points to a later paper. The multivalued nature of 
the dispersion surface is not simply a mathematical 
abstraction, but has been shown to be a physical 
reality in several laboratory experiments,—!^ and 
has been identified in several plasma descriptions, in¬ 
cluding gyrokinetics— In this appendix we identify 
two examples of exceptional points in the Vlasov- 
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Maxwell system, one at large MHD scales and a sec¬ 
ond at parallel wavelengths near the proton inertial 
length. 

The interchange of large scale fast and slow 
Vlasov-Maxwell waves )^i^^'^^ analogous to the inter¬ 
change between the MHD fast and slow modes, is 
due to an exceptional point at {/3p,0) (1.3,30°), 

where 6 is the angle between k and Bq. This excep¬ 
tional point, illustrated in the three panels of Fig¬ 
ure (HJ exists for large wavelengths with kpp < 1.0. 
To encircle the exceptional point, we follow the com¬ 
plex frequency solution for the fast and slow modes 
along four paths in (/3p, 9) space: from (I) (1.25,45°) 
to (1.58,45°), (II) (1.58,45°) to (1.58,15°), (III) 
(1.55,15°) to (1.25,15°), and (IV) (1.25,15°) back 
to (1.25,45°). All other plasma parameters for the 
system are fixed: kpp = 10“^, T±/T\\\p = Tj_/T|||e = 
1.0, Tp/Te = 1.0, and Vtp/c = lO”"^. The frequency 
and damping rates for the two modes along these 
four paths are presented as functions of either /3p 
or 9 in the top row of Figure HI We see that after 
following the prescribed path the two modes have 
exchanged complex frequency solutions. The con¬ 
tinuous nature of this mode exchange is shown in 
the lower left panel of Figure HI in which the same 
solutions for the fast and slow modes are plotted in 
complex frequency space over paths I-IV. This ex¬ 
ceptional point exists for wave lengths up to fepp ~ 1, 
as is shown in the bottom right panel, where the lin¬ 
ear solutions are calculated using the same paths in 
(/3p,0) with kpp = 0.1 (dashed lines) and 1.0 (dot- 
dash). For the latter, smaller scale case, the excep¬ 
tional point has disappeared and the fast and slow 
modes are no longer connected by a continuous vari¬ 
ation of j3p and 9. 

While the large scale (/3p, 9) exceptional point dis¬ 
appears at small wavelengths, a variety of other ex¬ 
ceptional points arise near k\\dp ~ 1 and kj_pp ~ 1 
connecting two, or more, modes through the varia¬ 
tion of other plasma parameters. An example of such 
a small scale exceptional point is found near /3\\p ~ 
9.0 and Tj_/T|| |p Ri 0.9. This exceptional point, illus¬ 
trated in FigureHl connects the Alfven and fast solu¬ 
tions. The values for the fast and Alfven solutions in 
complex frequency space along the paths connecting 
(/3||p,TV%) = (7.0,1.0), (7.0,0.7), (10.0,0.7), 
and (10.0,1.0) are calculated for a set of fc||dp G 
[0.2,1.0]. All other plasma parameters are fixed 
at k±pp = 10“^, T||p/T||e = 1, Tj_e/T\\e = 1 and 
Vt\\p/c = 10“"^. The value of k\\dp is indicated in 
Figure H] in color. Unlike the large scale (/3p,0) 
exceptional point, which persists over several or¬ 
ders of magnitude in wavelength k, this fast-Alfven 


exceptional point is narrowly restricted to k\idp G 
[0.6,0.9]; for scales above and below this range, the 
Alfven and fast solutions remain distinct. The ef¬ 
fects of this exceptional point can be seen in the 
left column of FigureHl where the /3||p = 10.0 and 
T±p/Ty < 1.0 Alfven solution with k\\dp = 0.6 
becomes unstable to the parallel firehose instabil¬ 
ity, the instability associated with the fast/whistler 
wave, a behavior not seen for any of the jdp < 10 
Alfven solutions. Note that the (dp = 10.0 fast solu¬ 
tion with fc||dp = 0.5, shown in the second column 
of Figure HI does not exhibit this interchange of so¬ 
lutions, illustrating the strong dependence of this 
exceptional point on k\\dp. 

Proper identification of linear solutions is com¬ 
plicated by the existence of exceptional points for 
Vlasov-Maxwell dispersion surfaces. The existence 
of these points allows for a continuous change in fre¬ 
quency and eigenfunction characteristics from one 
linear solution to another, allowing for possibility of 
coupling between distinct linear solutions and com¬ 
plicating any attempt to understand the nature of 
turbulent energy transfer and dissipation in terms of 
the nature of the underlying linear fluctuations. Un¬ 
derstanding the behavior of these exceptional points 
for collisionless plasmas and constraining their im¬ 
pact on the solar wind will be left to later work. 
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FIG. 5. Plots of uj/Q.p, 7 /fIp, 5'Q/5B, and 5n/no (in descending rows) for the four normal modes associated with 
proton temperature anisotropy instabilities (columns) as a function of T±p/T^^p. Values are plotted for two decades 
of (?||p, ranging logarithmically from 0.1 (blue) to 10.0 (red), with 1.0 given in black. Each of the modes have fixed 
values of fcxpp and fcydp. In the second row, stable (unstable) modes with 7 < 0 (> 0) are plotted with dashed (solid) 
lines. 
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FIG. 6. Plots of w/fc||t;A (top row), |(5i32,|/|(5i3y| (center), and {lo/ k\\VA)5hy/5Uy (bottom) for Alfven waves along 
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FIG. 7. Dispersion relations ij{k±pp)/Up for stable 
Alfven waves along the line of critical balance (dash- 
dotted lines) and for nnstable proton cyclotron and par¬ 
allel firehose modes with y/flp > 10“^, respectively illus¬ 
trated by the partially overlapping red and blue shaded 
regions. Plasma parameters are set to Tx^/Typ = 2.68 
(red), 2xp/T||p = 1.00 (black), and T±p/T^^p = 0.27 
(blue) with /3||p = 1.7. 
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FIG. 8. The transition between slow (blue lines) and fast (red) modes due to a large-scale exceptional point. The 
frequency ojjkvA and damping rate 'yjkvA for the two modes as a function of jSp and 6 along the path prescribed in 
the text are plotted in the top row. The parametrized path for these two modes in complex frequency space is given 
in the lower left panel for kpp — 10~^, and in the lower right panel for kpp — 0.1 and 1.0. 
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FIG. 9. Parametrized solutions in complex frequency 
space of the Alfven and fast modes along paths of /3||p 
G [7.0,10.0] and Txp/Typ G [0.7, 1.0], with the value of 
fc||cip G [0.2,1.0] given by the color bar. The mode con¬ 
version between the fast and Alfven modes for k\\dp G 
[0.6, 0.9] is indicative of the presence of a small-scale ex¬ 
ceptional point as described in the text. 
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